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Abstract 



Quantum theory of geometry, developed recently in the framework of non- 
perturbative quantum gravity, is used in an attempt to explain thermody- 
namics of Schwarzschild black holes on the basis of a microscopical (quantum) 
description of the system. We work with the formulation of thermodynamics 
in which the black hole is enclosed by a spherical surface B and a macroscopic 
state of the system is specified by two parameters: the area of the boundary 
surface and a quasilocal energy contained within. To derive thermodynamical 
properties of the system from its microscopies we use the standard statisti- 
cal mechanical method of Gibbs. Under a certain number of assumptions on 
the quantum behavior of the system, we find that its microscopic (quantum) 
states are described by states of quantum Chern-Simons theory defined by 
sets of points on B labelled with spins. The level of the Chern-Simons theory 
turns out to be proportional to the horizon area of black hole measured in 
Planck units. The statistical mechanical analysis turns out to be especially 
simple in the case when the entire interior of B is occupied by a black hole. 
We find in this case that the entropy contained within B, that is, the black 
hole entropy, is proportional to the horizon surface area. 
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I. INTRODUCTION 



The statistical mechanical method of Gibbs, developed in his celebrated Elementary 
Principles of Statistical Mechanics (1902), turned out to be the most efficient tool for ex- 
plaining the bulk properties of matter on the basis of its microscopies. So far, self-gravitating 
systems have resisted the application of this method. In this paper we apply the statistical 
mechanical method of Gibbs to such systems. The aim of the paper is to try to derive 
thermodynamics of gravitational systems using methods developed by Gibbs. 

Our discussion is restricted to a simple gravitational system whose thermodynamics is 
widely studied: a spherically symmetric uncharged black hole. Although this paper deals 
only with a special example of gravitational system, there is a natural way to generalize the 
ideas presented to other types of black holes (see [[[[]). 

As a basis of our discussion we take the thermodynamics of self-gravitating systems 
in the form it is formulated by Martinez ||. In this formulation the gravitational system 
composed of a spherically symmetric black hole is characterized by the surface area A = AirR 2 
of a two-dimensional spherical boundary surface B (located at r = R) that encloses the 
black hole, and a quasilocal energy E contained within. Following Martinez 0, we use the 
quasilocal energy of Brown and York |5| as E. In equilibrium states the system is completely 
characterized macroscopically when these two variables are specified (see 0). 

The expression of the entropy function S as a function S(A, E) of extensive variables is 
called fundamental equation of a thermodynamical system. Once known, the fundamental 
equation contains all thermodynamical information about the system. For the system of our 
interest the fundamental equation is given by 

S(A,E)=4nE*(l-^)\ (1) 

where 'radius' R is used as an extensive variable instead of the area A. Due to the spherical 
symmetry of the system one can use R and A interchangeably. 

The fundamental equation ([[]) is derived from Hawking's semiclassical expression ||] 
for the temperature of black hole radiation. It is our aim in this paper to try to explain 
the fundamental equation (|1|) on the basis of the microscopies of the system. To derive the 
fundamental equation ([I]) we apply the Gibbs' method. The method of Gibbs, when adapted 
to our case, is to construct a function Q(a, (3) called the statistical sum of the system, which 
is a function of intensive parameters of the system (here a stands for the product j3p of the 
inverse temperature (3 and the quantity p that plays a role of the 'surface pressure'; see 0). 
The statistical sum Q is given by 

Q(a,(3) = Tre- aA ~^, (2) 

where A and E are quantum mechanical operators corresponding to the classical quantities 
A, E. The statistical sum (|2|) contains all thermodynamical information about the system. 
In particular, the fundamental equation of the system can be obtained from Q(a,/3) by 
means of well-known thermodynamical relations. 

As a basis of microscopical description of our system we use the one given by non- 
perturbative quantum gravity based on the loop representation |J (see also |J). An im- 
portant modification arises, however, because our system has a boundary. In this paper we 
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propose a construction of the space of quantum states which takes into account the presence 
of the boundary. Having found quantum states of the system, to calculate the statistical 
sum Q one needs to construct quantum operators corresponding to the variables A, E, and 
find the corresponding spectra. The latter problem in the case of the area operator has been 
resolved successfully in ||. However, there exists only rather tentative results concerning the 
quasilocal energy operator. Surprisingly, even without knowledge of the operator E it turns 
out to be possible to analyze the thermodynamics of the system in a special case. Namely, 
we are able to complete the statistical mechanical analysis in the case when the boundary 
surface B that encloses the system coincides with the horizon surface of black hole. The 
analysis of this case constitutes the main result of the paper. 

We wish to note that our work uses some of the ideas presented in earlier works 0HH 
on black hole thermodynamics performed within the loop approach to quantum gravity. 
The important difference between the present and earlier works is that we propose to use 
Chern-Simons theory to describe states of the black hole quantum mechanically. 

The organization of this paper is as follows. Section describes the space of quantum 
states of our system. Section [TTT] constitutes the core of the paper. It contains a statistical 
mechanical analysis of the case when the whole interior of B is occupied by a black hole. 
We conclude with the discussion of the results obtained. 

II. QUANTUM STATES 

In this section we describe physical quantum states of our system. We give only main 
points of the construction. For details see [fflj . 

Let us recall that our system consists of black hole of a certain mass M enclosed by 
a spherical boundary surface B (located at r = R,R> 2M). Macroscopically the system 
is completely characterized by two variables, that is, the area A of B and the quasilocal 
energy of Brown and York E contained within, that are both defined on B. The quantities 
A, E completely define (via vacuum Einstein's equations) the state of the gravitational field 
within B. Therefore, one can think of a macroscopic (classical) state of our system as of a 
state of the gravitational field inside the surface B. 

We shall describe microscopic (quantum) states of the system as states of the quantized 
gravitational field within B. To construct these states we need to solve two problems, in 
complete analogy with the classical case. The first problem is to find the space of kinematical 
states; these are analogs of 'geometries' of the classical case. The second problem consists in 
imposing 'the quantum Einstein's equations'. Solving these equations one finds the so-called 
physical states. 

A major progress has been made on the first of these problems. Following the pioneering 
work m of Rovelli and Smolin, a mathematically well-defined theory of 'quantum geometry' 
has been constructed (see |],[| for the most recent account). However, a general treatment 
of 'the quantum Einstein equations' is still a matter of much controversy, despite the recent 
progress on this front fp~Q| ] . 

In this paper we propose a solution to the second problem for our special case, and 
find the physical states. To describe the main idea of our treatment we need a simple 
observation from the classical theory. Let us note that in the classical case the geometry 
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within B is completely determined (via vacuum Einstein's equations) once the geometry on 
the surface B is specified. Indeed, for our special case of static Schwarzschild spacetime 
dynamics (in the sense of the Hamiltonian formulation) is trivial. Therefore, the problem 
of solution of Einstein's equations reduces to the problem of solution of constraints, which 
is an elliptic problem. This means that, to find the geometry within B, one should specify 
an (allowable) geometry on B (or, so to say, boundary conditions). This determines (via 
constraint equations) the geometry within the entire B. 

Following the spirit of the correspondence principle advocated by N.Bohr one can expect 
that, whatever precise form of equations of quantum dynamics is, these equations will keep 
some key properties of their classical analogs. For example, it is natural to expect that 
solutions to quantum dynamical equations describing static (in some appropriate sense) 
geometry within B should be completely determined by a state of geometry on the boundary 
B, as this happens in the classical case. We assume that (still unknown) equations of 
quantum dynamics possess this property. The assumption implies that physical states of 
our system are labeled by quantum states of geometry on the boundary surface B. Note 
that our assumption does not tell us anything about non-stationary cases for which the 
corresponding classical property does not hold. 

Thus, the analogy with the classical case tells us that equations of quantum dynamics 
should be such that, once a state of quantized gravitational field on B is specified, one finds 
a state of quantized gravitational field within the entire B solving these equations. However, 
even in the classical case there does not exist a solution for any choice of 'boundary data'. 
We expect, therefore, that in the quantum case not all states of quantum geometry on the 
surface B will give rise to physical states of the system. Moreover, in the classical case we 
need to know only certain minimal set of data to find the geometry within B. We expect to 
have an analog of this in the quantum case. 

Thus, to construct quantum states of the Schwarzschild black hole let us discuss the 
classical case in more details, and find a minimal set of surface data that completely specifies 
a macrostate of the system. We describe the gravitational field using Ashtekar variables 
[l]J. The geometry of spatial hypersurfaces S 4 is described in this formalism by a pair 
of (canonically conjugated) variables E a , A a . Let us now introduce fields describing the 
geometry on B. First, we introduce the pullback a of the connection A on B. Next, let us 
denote by e the two-form that is the pullback on B of the dual to the densitiezed triad E a 
two- form e a b c E°. 

On solutions of Einstein's equations the surface fields a, e are not independent. In the 
Schwarzschild spacetime there exists a simple relation between the curvature two-form of 
the connection a and the two-form e. Taking from [O] explicit expressions for a, e in the 
Schwarzschild spacetime, one can find that 

2M 1 
f ~ R R 2 

Note that throughout this paper we use the units in which h — G — c — 1. 

We take (|3]) as the sought surface constraint equation. Thus, we view two fields e, a on 
the surface B, where the connection a is such that its curvature two-form is related to e via 
/ = if)e, ip being a number, as specifying a macroscopic state of our system. Indeed, it is 
easy to see that classically this set of data carries all information necessary to reconstruct 
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the geometry within the entire B. Let us note that the two form e carries information about 
areas of regions on B, and, in particular, the area of the entire B can be determined once e 
is known. If, moreover, one is given that the curvature two-form / of the connection on B is 
related to e via / = ipe, one can find the mass of the corresponding Schwarzschild solution 
from ([|). Thus, a set of fields e, a together with the relation (|3]) between them carries all 
information about a macrostate. 

Following our method of analogy with the classical case, we now simply have to find how 
to describe states of fields e, a quantum mechanically, and how to impose the constraint (|3]) 
in the quantum case. The states we find this way will label the physical states of our system. 

We are lucky, for there exists a complete description of states we look for. Quantum 
theory of geometry |J provides us with a complete description of states of quantized e 
field. Conformal field theory JR] tells us what are quantum states of connection subject 



to a quantum analog of constraint (y) on a two-dimensional surface. The corresponding 
quantum states, although in a different context, have already been used in non-perturbative 



quantum gravity by Smolin fL5| . In this paper we simply give a description of these quantum 
states in the amount we need for our statistical mechanical treatment. 

The states we are going to describe are invariant under diffeomorphism transformations 
that are tangent to the boundary surface B, gauge invariant, and satisfy, in some precise 
sense, the quantum analog of the constraint We give a description of a basis in the space 
of quantum states that is formed by eigenstates of operators measuring areas of regions on 
B. This basis is especially convenient for our purposes. A basis quantum state is labeled by 
a set of points on the surface B, which we, following ||, shall call vertices. It is important 
that two states that differ only in position of vertices on B should be considered as a single 
state if one deals with diffeomorphism invariant states. Each vertex v is labeled by a set 
of quantum numbers. These include: (i) spin (half-integer) that enters the interior of B 
starting at the vertex v; (ii) spin j" that enters the surface B on an edge lying up the surface; 
(iii) spin j* that stays on the surface (see Fig. |IJ). Because of gauge invariance quantum 
number j* can take values only from the range |j" — j$\ < il < j% + jy- Note also, that 
because of gauge invariance not all sets of spins are allowed ||. Namely, in the case of a 
closed surface, which we consider here, the half-integers should satisfy the additional 

condition that the sums J2 v ^b 3v anc ^ ^veBJv over au vertices in B are integers. 

The set of quantum numbers we have described can be thought of as specifying a state 
of e field on B. As it is described by conformal field theory, for each choice of these quantum 
numbers there exists a finite number of states of the quantized connection a on B that satisfy 
the constraint ([3D. These states are the states of quantum Chern-Simons theories defined 
by a set of spins {j^}- Different possible states of the quantized connection a are labeled by 
additional quantum numbers. 



It turns out (see, for instance, |16| for a discussion of this point) that the zero vector is 



the only vector in the physical Hilbert space of states of quantized connection a, unless the 
representations of SU(2) labeled by spins j* all satisfy a certain condition of integrability. 
It is shown, for example, in that for the case when the constraint equation is 



the integrability condition is simply that all spins j* satisfy the inequality j* < k. Here k is 
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required to be an integer. Comparing ([|) with the relation (J3|) we find that in our case 



k ~ A j2M/Wy (5) 

It is known that representations satisfying the above integrability condition are, in fact, 
representations of the quantum group SU(2) q . The number k is known as level, or coupling 
constant of the corresponding Chern-Simons theory ||15|| . It is related in a simple way to the 
deformation parameter q (see, for example, Hl5|). 



It is not hard to see, that the level k takes the minimal (for a fixed area A) value A 
when the surface B coincides with the horizon surface of black hole. When the mass of black 
hole inside B decreases, the level increases, as it is obvious from (|5|). Thus, the level k is 
proportional to the area of the surface B (measured in Planckian units) when this surface 
coincides with the horizon surface of black hole, and this value is the minimal possible value 
once A is fixed. This has an important consequence to us, for this means that to large (in 
Planckian units) black holes correspond large levels k. It is well-known that in the limit 
of large k the dimensions of the quantum group representations go into dimensions of their 
classical analogs. Also, there exists a simple formula for the number of different states of 
quantized connection in the limit of large k. Roughly speaking, a basis in the space of states 
of quantum connection a is labeled by different ways that the spins j* can be combined 
consistently according to the rules of addition of angular momentum of the quantum group 
SU(2) q . In the limit of large k this number is given by (see |15| ) 

nc# + !)• ( 6 ) 

V 

One can expect that statistical mechanical methods can only be applied to macroscopical 
objects that are composed of many 'elementary' excitations. For our geometrical system this 
means that statistical mechanical description is presumably legitimate only when applied 
to systems large as compared with Planckian scales. As we have seen, for large black holes 
the level k of the quantum theory is large, and one can use the simple formulas that hold 
in the limit of a large k. Thus, we shall not keep track of the fact that k is finite. This can 
be shown to be legitimate, for the corrections that appear when one takes the finiteness of 
k are negligible. 
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FIGURES 




FIG. 1. A typical quantum state of our system is labeled by a set of vertices, with quantum 
numbers (spins) attached. The spins j should satisfy the condition j < k. For every set of vertices 
on B (and quantum numbers) there exists only a finite number of states of quantum connection, 
which can be thought of as the number of different ways in which the spins j\ can be combined 
consistently according to the rules of addition of angular momentum of the quantum group SU{2) q . 

The states we have described (see Fig. |] for a picture of a typical quantum state) are 
eigenstates of operators that measure areas of regions on B. In particular, the eigenvalues 
of operator An that measures the area of a region 1Z on B are given by p| 



A n (S) = 167T7 5J « V 2 #(# + !) + WW + + !)> 



(7) 



where we have introduced the notation S for the basis quantum states. The sum here is 
taken over all vertices of S that lie in the region TZ. We have taken into account in (0) the 
fact that states S are gauge invariant. The quantity 7 in ([/j) is the parameter that appears 



in the loop quantization of general relativity, as it is discussed in ||17|| . In the case of closed 
surface, which we consider here, not all of the eigenvalues (^) are eigenvalues of the operator 
A that measures the area of the entire B. Namely, in the case when B is closed gauge 
invariant quantum states are those that satisfy the additional condition that we mentioned 
above. 

Unfortunately, there does not exist yet such a complete description of operator E that 
measures the quasilocal energy. However, as we shall see in the next section, a partial 
statistical mechanical analysis is possible even without knowledge of E. 



III. STATISTICAL MECHANICAL ANALYSIS 

The statistical mechanical method of Gibbs consists in constructing the statistical sum 
over all different quantum states of the system. Since a macroscopic state is specified by the 
two extensive parameters: the area A and the quasilocal energy E, we have to introduce two 
conjugate intensive parameters that we shall denote by a, (3 correspondingly. The statistical 
sum is a function of the intensive parameters 

Q(a,f3) = Tre- aA -^. (8) 
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The statistical sum can not be calculated, unless the operator E is known. Let us note, 
however, that for two values of the parameter (3, namely, for (3 — 0, oo, the calculation is 
possible. Indeed, for (3 = a precise form of the operator E simply does not matter. In the 
latter case only the states of zero energy survive in (^). 

It is obvious that the latter case describes the zero energy state of the system (it cor- 
responds to the zero temperature, that, as we know, is proportional to the inverse (3). In 
other words, in the case (3 = oo our system is simply a region of flat spacetime enclosed 
by the surface B. When one decreases (3 (increases the temperature) the mean value of 
energy E will increase. One finds that the value (3 = (infinite temperature) corresponds 
to the largest possible value of energy for fixed A. But, we know that the maximal (for 
a fixed area A) possible 'amount' of quasilocal energy is contained within B when all the 
interior of B is occupied by a black hole. Thus, the case of (3 = will describe the situ- 
ation when B coincides with the horizon surface of black hole. This is in agreement with 
the thermodynamics described in 0. Indeed, (3 becomes zero when the surface B coincides 
with the horizon surface of black hole (see 0]). Thus, we find that some thermodynamical 
functions describing black hole can be calculates even without knowledge of E. To calculate 
thermodynamical functions in this case, we simply have to put (3 = in (||). 

To find Q(a,P) in the case (3 = let us use the basis of eigenstates of A described in 
the previous section. It is convenient to think about vertices on the surface B as about 
imaginary particles. Different particles correspond to different values of quantum numbers 
j u , j d , j t . There can be any number of particles of each sort in a state. It is important to take 
into account the fact that, according to formula (|B|), each particle 'carries' the degeneracy 
(2j* + 1) (that is, the addition to a state of a particle increases the dimension of space of 
states of quantized connection). The sum over all states becomes the sum over numbers of 
particles 

«(<•) - II £ mr exp <-„ n r 4(F)) = FJ j _ ^(-^(r)) ' < 9 > 

where we have denoted by T a sort of particles (r simply labels a set of spins j u ,j d ,j t ), 
by d(T) := (2j* + 1) the degeneracy carried by a particle of sort T, and by nr the number 
of particles of sort T in a state. The summations over nr run from zero to infinity, A(T) 
denotes the contribution to the area of B from a particle of sort T. 

For simplicity, we have not taken into account in (|^) the fact that not all eigenvalues 
(0) are eigenvalues of area operator for a closed surface. In our language of particles this 
means that not all combinations of particles can be realized. Thus, strictly speaking, the 
summations over numbers nr of particles in (0) are not independent. Let us forget about 
this for a moment and proceed with our analysis. We shall discuss the consequences of this 
restriction later on. 

We have found that the state of the system when the entire interior of B is occupied by 
a black hole corresponds to (3 = 0. Thus, such a state of the system is described by a single 
thermodynamical parameter a. The mean values of the area A and the energy E, that can 
be expressed as derivatives of In Q(a, (3) at (3 = 0, become functions of this single parameter. 
One could exclude a, and express all thermodynamical functions in terms of, for example, 
A. In other words, a macroscopic state of a Schwarzschild black hole is specified by a single 
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parameter. We are particularly interested to find the entropy function of the system as a 
function of A. Clearly, it may serve as a first test to our method whether it predicts the 
Bekenstein-Hawking formula for the entropy of black hole. Namely, the thermodynamics 
described in tells us that, when the whole of B is occupied by a black hole, the entropy 
'contained' within B is simply the Bekenstein-Hawking entropy of black hole 

Sbh = \a (10) 

To find the entropy as a function of A one should, in principle, find S(a) and A(a). 
Then, one would exclude a and find S as a function of A. This problem can, in principle, 
be solved; however, there exists a simpler way to get the dependence S(A). We note that 
we expect the Bekenstein-Hawking formula ( |10|) to hold only for large (as compared with 
Planckian area) areas A. Indeed, from the point of view of quantum geometry only large as 
compared with the Planckian scale black holes can be considered macroscopical, and thus, 
the thermodynamical description can be applied only to such black holes. Thus, we are 
interested in the dependence S(A) predicted by the statistical mechanics only for large values 
of A. One can change the mean value of area changing the intensive parameter a. Some 
values of a correspond to large mean values of A, and we are to find the thermodynamical 
functions for these a. 

From the formula 

dlnQ(a) 

da 

we find for the mean value A(a) of the area 

A(a) = 5>(r)/(r), (12) 

r 

where we have introduced the function 

= d(r)exp(-aA(r)) 
J{ ' l-d(r)exp(-aA(r))' 1 ' 

It is easy to see from ( O ) that f(T) plays the role of the mean number of particles in the 
'state' T for fixed value of parameter a. Large values of A(a) correspond to a case when at 
least for some T the function f(T) takes large values. 

It can happen that the second term in the denominator of ( |13|) for some T is close to 
unity. This will correspond to a large number of particles of the sort T, and, therefore, to a 
large mean area. Let us see to what values of a this corresponds. 

As it was noted in |7|], there exists only few possibilities. First, it can happen that the 
thermodynamical functions such as Q(a), A(a), S(a) get large values only when a — > 0. It 
follows then from the Euler relation that in the thermodynamical limit (large areas) the 
dependence S(A) of the entropy on the area is different from the linear one (in fact, S(A) 
grows slower than A). The second case is when the thermodynamical functions diverge as 
a goes to some 'critical' value a cr . Then, as it is discussed in J7|, the Euler relation implies 
the linear dependence S(A) = a cr A. Finally, there exists the possibility that the statistical 
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sum Q diverges for all a, which means that the density of states of the system grows faster 
than exp A. 

It is not hard to show, that in our case the second possibility is realized, and, therefore, 
the dependence S(A) of the entropy on the area (for (3 = 0) is linear. Indeed, let us consider 
the function d(F) exp — a A(F) as a function of T and a. Let us recall that T stands for 
three quantum numbers j u ,j d ,j t - It turns out to be convenient to use instead of j u ,j d ,j t 
the following linear combinations of them 

I:=f + j d , 

J--=f-J d , (14) 
K ■= j u + j d — j u+d 

Then, as it is not hard to see, /, J take all positive integer and half-integer values, and 
K runs from zero to 2min(j u , j d ) (K also takes both integer and half-integer values). The 
contribution to the area from a particle of sort F written as a function of this new quantum 
numbers becomes 

A(F) = A(I, J, K) = 87T7V J 2 + 1 + 21 K + K - K 2 (15) 

The degeneracy that is 'carried' by a particle of sort T becomes 2(1 — K) + 1. As a straight- 
forward analysis shows, for all values of a < a cr = 1.138/87T7 t ne second term in the denom- 
inator of (fL3|) is less than unity for all types of particles T (for all non-zero values of quantum 
numbers I,J,K). When a ~ a cr , the maximal value of the function d(F) exp — a A(T) in 
the allowable range of the quantum numbers /, J, K is just the unity. This is realized for 
the following values of quantum numbers: j u = j d ;j u + j d = j l = 2. 

Thus, what we find is that for the value of a close to a cr the number of particles f(T) in 
one of the 'states' F = F cr becomes large. This means that the thermodynamical functions 
diverge for a = a cr . As we have said above, this implies that for large values of A the 
dependence of the entropy on the area is linear 

S = a cr A, 

a cr = 1.138/87T7 (16) 

Note that we have obtained this result analyzing the behavior of f(F) as a function of F 
and a. It is crucial for the result that our particles 'carry' the degeneracy d(F), for it is 
this degeneracy in ([13]) that makes the number of particles in one of the states diverge for 
a finite (non-zero) value of a. 

As we have said, when a approaches a cr the number of particles in one of the states 
T = r cr goes to infinity. It is interesting to compare this with the phenomenon of Bose- 
Einstein condensation. In the case of an ideal Bose gas a macroscopic fraction of particles 
that compose the gas end up on the ground energy level when the temperature of the gas 
reaches certain critical temperature. Let us note that the number of particles in the case of 
a Bose gas is fixed. In our case particles are imaginary creatures, and their number is not 
fixed. As we have found, when the intensive parameter a approaches the critical value a cr , 
much larger number of particles end up in the state F cr than in all other states. Thus, what 
we find is very similar to Bose-Einstein condensation; the important difference is, however, 
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that the number of particles in our 'gas' is not fixed. Indeed, the fact that the number of 
particles is not fixed can be shown to be crucial for our result of linear dependence of the 
entropy on the area to hold.Q 

When a is close to a cr only a negligible fraction of particles is in states different from 
r cr . Thus, one can write 

A = A(T cr )N, (17) 

where we have introduced the number N of particles in the state r cr . 

Let us now recall that all results of this section were obtained neglecting the fact that, 
in the case of a closed surface B, not all sets of spins on B correspond to physical states. As 
we have mentioned above, there exists a simple condition on the total sums of spins that 
enter and leave the surface. The way to take into account this condition is to exclude from 
the statistical sum all terms with spins not satisfying the conditions. Generally, such an 
operation may significantly change the behavior of all thermo dynamical functions that are 
derived from the statistical sum. It is not hard to see, however, that our condition is very 
weak in the sense that the results we have obtained for the case of an open surface continue 
to hold when one considers closed surfaces B. Here we present only a heuristic argument in 
the support of this, leaving a more rigorous treatment for another occasion. 

It is not hard to see that the results fljif) , (|T7|) stem from one and the same fact that 



we have discovered about the microscopic state describing black hole. We saw that a state 
describing a large black hole is very special in the sense that a macroscopic fraction of 
'particles' resides in one and the same 'state'. Then the total area is given simply by the 
number of particles times the contribution from an individual particle, as it is described by 
(|T7j). The entropy, on the other hand, is given by the logarithm of the number of states of 
quantized connection 

S = In d{T cr ) N = N In d{T cr ). (18) 
Combining this equation with (|T7|) we find that 



A(T cr ) 8n-fV2 1 ; 

Since j u ,j d in this state happen to have one and the same spin 1, the conditions we have 
mentioned above are satisfied in any state describing a large black hole. Thus, in the 
thermodynamical limit (large black holes) the results we have obtained continue to hold 
even when one imposes the conditions arising in the case of a closed B. 



A similar phenomenon, when the entropy of the system grows linearly as a function of energy, is 
known to occur in ordinary thermodynamics for systems for which the number of particles is not 
fixed. The author is grateful to P. Aichelburg from whom he learned about this fact. 
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IV. DISCUSSION 



We have found that the microscopic states of Schwarzschild black hole can be described 
by states of SU(2) Chern-Simons theory, which are defined by choices of vertices and spins 
on B. We have found that the integer level of Chern-Simons theory is proportional to the 
horizon area of black hole (see (^j)). Thus, large black holes correspond to large levels k. We 
have used this description as the basis of our statistical mechanical analysis. 

Although we have not reached our goal of explaining the fundamental equation ([!]) on 
the basis of quantum microscopies of the system, we were able to complete the statistical 
mechanical analysis of the case when the entire interior of the system is occupied by a 
black hole. We have found that in this case the entropy contained within B is proportional 
to the area of the boundary B, with the proportionality coefficient given by The 
proportionality coefficient between black hole entropy S and the horizon area A turns out 
to be a function of the parameter 7. 



As it is discussed in |[L7|| , the parameter 7 is a free parameter that arises in the loop 
quantization of general relativity. Unless a value of the parameter 7 is fixed by some in- 
dependent considerations, a comparison of the dependence S(A) predicted by our analysis 
with the Bekenstein-Hawking formula (|TC|) is not possible. Thus, the predicted dependence 
([L6|) by itself does not provide us with a test of our approach. Note, however, that the 
approach presented gives us not just the dependence fll6|) . The considerations of section |T| 
led us to the conclusion that the statistical state of a black hole should be described by a 
density matrix 

b = 1 e~ aA 

where Q(a) is the statistical sum (Q). The entropy was then defined as S = Tr(plnp). 
However, with the density matrix p in one's disposal one knows much more than simply the 
entropy. For example, the knowledge of the density matrix allows one to analyze properties 
of the black hole radiation spectrum fl8| . Results of such an analysis, together with the 



result (|T6|), 

can serve as a test of the validity of the approach presented. 
The other result of this paper is the discovery of the fact that the state of the system 
that corresponds to a macroscopically large black hole is realized for values of the intensive 
parameter a that are close to the critical point a cr . This is quite similar to what one finds, 



for example, in the sum over lattices approach to quantum gravity in two dimensions |T9 
One finds that, in order to go into a macroscopical limit, one needs to tune values of the 
parameters to a critical point. It is interesting that in our case, when the parameter a goes 
to a critical point, the system undergoes a phase transition similar to the phenomenon of 
Bose-Einstein condensation. We find that a macroscopical fraction of elementary excitations 
of geometry is in one and the same quantum state when we are dealing with large black 
holes. The state where most of the particles 'condense' to happen to have the following 
quantum numbers: j u = j d = 1. 

Let us conclude summarizing the open problems of our approach. One of the most 
important of such problems is to develop a theory of quantum dynamics of geometry, and to 
answer the question whether the assumptions about dynamics made in this paper are true. 
This is, however, a difficult problem that may require mutual efforts of a large number of 
researchers in the field. 
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Probably a simpler, but not less important problem is to construct the operator E 
corresponding to the quasilocal energy. This would allow one to analyze the case of an 
arbitrary mass contained within B and find the fundamental equation of the system. Work 
is in progress in this direction. 

Another important open problem, related to the problem of construction of the energy 
operator, is to understand how the different mathematical techniques used in conformal 
field theory and in non-perturbative quantum gravity can be reconciled in a mathematically 
rigorous construction of quantum states. Work is in progress also in this direction Q. 
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